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ABSTRACT 



Controversy regarding the correctness of a test for alias- 
ing proposed by Hinich and Wolinsky [5[ has been surpris- 
ingly long-lived. Two factors have prolonged this contro- 
versy. One factor is the presence of deep-seated intuitions 
that such a test is fundamentally incoherent. Perhaps the 
most compelling objection is that, given a set of discrete- 
time samples, one can construct an unaliased continuous- 
time series which exactly fits those samples. Therefore, the 
samples alone can not show that the original time series was 
aliased. The second factor prolonging the debate has been 
an inability of its proponents to unseat those objections. In 
fact, as is shown here, all objections can be met and the 
test as stated is correct. In particular, the role of station- 
arity as knowledge in addition to the sample values turns 
out to be crucial. Under certain conditions, including those 
addressed by the bispectral aliasing test, the continuous- 
time signals reconstructed from aliased samples are non- 
stationary. Therefore detecting aliasing in (at least some) 
stationary continuous-time processes both makes sense and 
can be done. The merits of the bispectral test for practical 
use are briefly addressed, but our primary concern here is its 
theoretical soundness. 

1. THE BISPECTRAL ALIASING TEST 

The domain of the discrete-time bispectrum is the two di- 
mensional bifrequency {wi,W2} plane. Assuming a real- 
valued discrete time series, the usual replication phenomenon 
dictates that all non-redundant information is confined to the 
square < u>i, u>2 < 7r. When one fully accounts for sym- 
metries, the non-redundant information in the bispectrum is 
confined to a particular triangle inside this square |[y, g]. 

This triangle naturally divides into two pieces. One piece 
is an isosceles triangle and is unproblematic. The other 
piece, somewhat unusual in shape, is the source of the con- 
troversy under discussion. Naive consideration of this trian- 
gle shows that it involves frequencies higher than the Nyquist 
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frequency and therefore must have something to do with 
aliasing. Hinich and Wolinsky considered this more care- 
fully and showed that the naive intuition is correct: if the 
discrete time series arises from sampling a stationary, band- 
limited, continuous-time process, and if the sampling rate 
is sufficiently rapid to avoid aliasing, then the discrete bis- 
pectrum is non-zero only in the isosceles triangular subset 
of the fundamental domain. Conversely, if the bispectrum 
of a sampled stationary continuous-time process is non-zero 
in the outer triangle, then the sampling rate was too slow to 
avoid aliasing. 

It should be clearly understood that there is no assertion 
that aliasing in general can be detected. The statement is 
not "if a signal is aliased, then the outer triangle will have a 
non-zero bispectrum." Rather, the assertion is the converse, 
"if the outer-triangle shows a non-zero bispectrum, the (un- 
derlying) continuous-time signal must have been aliased." 

At one level, this result is obvious and, in fact, the result 
was initially so-regarded [ph- However, doubt soon arose. 
Perhaps the most important source for suspicion is the argu- 
ment based on reconstruction alluded to above. 

In light of this objection, one is led to reconsider the 
association of the outer triangle with aliasing. One can take 
the position that there is no relation, as in [ph. One can 
decide that something is aliased, but that it is the bispectral 
estimator rather than the signal. There is some plausibility 
to this claim, for the frequencies that are involved in the 
outer triangle are u>x, u>2, and u>\ + uj-2 — 27r. This seems to 
be the position of Pflug et al. [H]. 

Or, one can try to delineate the conditions, if any, under 
which the test makes sense. This was done by Hinich and 
Messer in 1995[Q]. They confirmed the validity of the orig- 
inal argument and stated its conclusions more carefully. In 
particular they conclude that a non-zero bispectrum in the 
outer triangle indicates a non-random signal or one of the 
following: 

• a random, but non-stationary signal ; 

• a random, stationary, but aliased signal, or; 

• a random, stationary, properly-sampled signal which 
violates the mixing condition. 



We believe that the analysis of Hinich and Messer, while 
entirely correct, did little to persuade the detractors of the 
test. In particular their analysis did not address the recon- 
struction objection and may have left the impression that the 
circumstances for which the test applies are unlikely to be 
met in practice. 

In this paper, we show that the reconstruction objection 
is far from fatal. We further establish that stationarity is the 
only property which is crucial to the test. Since this prop- 
erty is required in order to define the bispectrum, one can 
legitimately apply the aliasing test whenever one is entitled 
to compute a bispectrum. Therefore the bispectral aliasing 
test is as theoretically sound as the bispectrum itself. 

2. THE SELECTION RULE AND BRILLINGER'S 
FORMULA 

The bispectrum, defined to be the triple Fourier transform 
of the third-order autocorrelation, reduces to a function of 
two frequencies since stationarity confines the spectrum to 
the plane through the origin of the frequency domain per- 
pendicular to the vector (1,1,1). 

di23(cs{ti,t2,t 3 )) = b(0Jl,U! 2 )S(uJl +uj 2 +w 3 ) (1) 

Another way of computing the bispectrum is to switch 
the order in which one does the Fourier transforming and 
the ensemble averaging. This leads to the following result. 

b(cji,uj 2 ) = {X(u>i)X(u)2)X(uJ3 = loi +lu 2 )) (2) 

If the process is bandlimited and X(io) — for |tt>| > it, 
then the bispectrum is confined to the intersection of the 
(1, 1, 1) plane and the 7r-cube (i.e. (o>i, u 2 ,u^) 6 (— 7r, 7r)<g> 
(— 7r, 7r) (g> (— 7r, 7t) ). The plane and its projection onto the 
(uji,u> 2 ) plane is shown in Figure 1. Upon sampling with 
unit time step, one obtains the usual replication in three di- 
mensions. (Doing everything in 3-dimensions and project- 
ing at the end keeps things simpler and makes it easier to 
avoid errors.) In particular, one gets that if the process is 
sampled at a frequency greater than twice the highest fre- 
quency component, then the bispectrum is confined to the 
replications of the tilted hexagon shown. 

The replication gives the discrete-time bispectrum b c V- 
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b(u>i,u) 2 ,oj 3 ). 



(3) 



where u>i = \ + 2?rk for integer k. 

Since the replication does not cause any overlaps, the 
outer triangle remains empty. This is the Hinich and Wolin- 
sky aliasing theorem. (Note that the outer triangle is equiv- 
alent to the bigger triangle with vertices (0, 7r, 0), (0, 0, it) 
and (0, 7r, 71") by symmetries. See [nl Bfl for details.) 




Figure 1 : The origin of the bispectral fundamental domain. 

3. THE RECONSTRUCTION OBJECTION 

Suppose we have a stationary process x(t) and we under- 
sample it by sampling at t € Z, Then by convolving x(t) 
with the appropriate sine function we get a reconstructed 
process x r (t). This new process will have exactly the same 
samples as the original process and therefore exactly the 
same sampled bispectrum: yet it is not aliased. Therefore 
for any process that is undersampled, we have another pro- 
cess producing an identical sampled process which is not 
undersampled, showing that that one could not possibly de- 
tect aliasing via the bispectrum computed from samples! 

The rub here is the fact that the reconstructed signal will 
not necessarily be stationary. Processes reconstructed from 
aliased samples of continuous-time signals are generally cy- 
clostationary but not stationary. Some aliased processes do, 
in fact, reconstruct into stationary processes. But in the class 
of stationary signals for which the bispectral aliasing test 
gives positive results, reconstruction from aliased samples 
produces non-stationary processes. 

To carefully illustrate this we will consider several sta- 
tionary processes generated by taking a periodic signal with 
period T and giving it a random shift 9 £ [0,T). First con- 
sider a simple cosine process, 



x(t) = cos(airt + 2ir9/T), 



(4) 



where a = 1.5, T = 4/3, and 9 is randomly chosen from 
[0,4/3). Upon sampling and reconstruction we get the co- 
sine process given by 



x r (t) = cos(a7rf + 2tt9/T), 



(5) 



where a — —0.5. The key idea is that the signal appears 
at the lower frequency as dictated by its replication into the 
fundamental region of the frequency space, (-7T, 7r), but its 
reconstructed phase is the same as the "source" component 
phase. 

Now consider 

x(t) = cos(airt + 2ira9) + cos(/3irt + 2ir/39), (6) 



where a = 1.0, (3 — 3.0, and 9 is chosen randomly from 
the interval [0, a -1 ). The reconstructed process one gets is 



x(t) = cos(a7rt + 2-KaO) + cos{{3irt + 27T/30), (7) 

(where a and /3 are the aliased frequencies). Although the 
phase terms 2na9 and 2ir(39 are still uniformly distributed 
over 2tt they now correspond to different time shifts so we 
no longer have a single shifted waveform. This process can 
easily be seen to be cyclostationary, but not stationary. See 
Figure 2 [g], Computation of the required expectations re- 
quires that one can "average over the ensemble." Since this 
stationary process is not ergodic, one can not get the result 
from a single realization of the process. It is at this point that 
some differences in perspective arise. Strictly speaking, in 
order to compute a bispectrum one must perform the ensem- 
ble average. A single realization does not suffice unless the 
process is ergodic. 




Figure 2: Cyclostationarity of signal reconstructed from 
aliased samples. The lower curve is the process envelope. 
The sampling interval is e. The upper curve is the sixth 
moment, chosen for ease of display. 

Finally, consider the process defined by 



x(t) = cos((10/20)7rf + 5 • 2ir9) -i 
cos((12/20)7rf + 6 • 2tt9) -+ 
cos((22/2O)7rt+ll-27r0) 



(8) 



where 9 is chosen randomly from [0,1). 

Because the phases of these components are in a fixed 
relation, this process has a spike in the bispectrum at wi = 
10/207T, uj 2 = 12/20-7T, i.e., in the outer triangle. The re- 
constructed signal is given by 



x(t) = cos((10/20)7rf + 5 • 2it9) -\ 
cos((12/20)7rf + 6 • 2x0) 4 
cos(-(18/20)7rt+ 11-2x6 



(9) 



which is not stationary. Therefore we have a signal with 
nonempty outer triangle whose reconstruction is not station- 
ary. This situation is exactly what the bispectral test implies 
happens whenever the outer triangle is nonempty. The loss 
of stationarity causes the Fourier transform of the triple au- 
tocorrelation to "move off" of the (1,1,1) plane. 

Therefore, if one knows (or is willing to assume) that 
the process which generated the observed samples was sta- 
tionary, one can rule out the unaliased reconstruction as the 
source of the samples. In a sense, the continuous time sig- 
nal reconstructed from aliased samples of an original time 
series is a "measure zero" object. This result is very surpris- 
ing to most people's intuitions. It is studied further in Vixie, 
Sigeti and Wolinsky ||. 

4. THE REPLICATION OBJECTION 

Upon looking at Equation one may observe that even if 
X(oj) = for |w| > 7T, sampling effectively fills in the spec- 
trum at higher frequencies. This is the basis for the objec- 
tion appearing in Swami \un. This concern is addressed as 
follows. While the spectrum does indeed fill out upon sam- 
pling, the undesired expectations remain zero. Consider a 
(statistically stationary) ensemble constructed by uniformly 
translating a periodic or finite-duration waveform x(t). Two 
operations are necessary to produce the discrete-time en- 
semble; a uniform shift in time over a period T, which in- 
troduces linear phase factors, and sampling, which produces 
spectrum replication. These operations do not commute: 
i.e., one wants to time-shift the waveform first and then sam- 
ple, rather than to shift its samples. For the shifted samples 
x s (t + 9) one finds 



d(x s (t + 9)) = e^d(x s (t)) 



(10) 



But for the sampled shifted waveforms x(t + 9)\ s , i.e., 
the waveforms needed to construct a stationary ensemble, 
the phase of the original signal is propagated to higher fre- 
quencies periodically rather than linearly. This difference 
leads to the vanishing of unwanted expectations. 

For example, consider the process given by the randomly 
shifted sum of unit amplitude cosine waves with frequencies 
at ri/20 (rad/s) where n takes integer values from 1 to 19. 
The sampled spectrum has components at u>\ = 107r/20, 
uj2 = 117I-/20 and lo^ = — 2l7r/20 but the average 



reduces to 
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(11) 



(12) 



where 9 is chosen with uniform probability from [0,1). This 
average vanishes. Therefore, the potential contribution in 
the outer triangle is zero because averaging kills it. This is 
in contrast to the case where the average is zero because the 
spectral amplitudes are themselves zero (as in the proof of 
the aliasing test). 

5. EMPIRICAL COUNTER-EXAMPLES 

Other objections to the test have been made. Frequently 
these objections involve a (purported) counter-example to 
the bispectral aliasing test. A particularly clear example is 
provided by Frazer, Reilly and Boashash [Q|. Here the au- 
thors do two things. They present an example of an aliased 
signal which the aliasing test fails to mark as aliased. The 
example is unproblematic: neither the aliasing test nor any 
aliasing test we are aware of will detect all aliased signals. 
It is not, however, a counter-example to the test. Since there 
is nothing in the outer triangle, the bispectral aliasing test 
makes no assertion regarding the presence of aliasing. 

The other example the authors provide is more interest- 
ing. It consists of a signal involving coupled sinusoids at 
Wi = 0.3125-Hz, u 2 = 0.25Hz and w 3 = A375Hz and 
the authors show that there is a peak in the outer triangle un- 
der conditions which rule out aliasing. As the authors note 
these frequencies sum to 1 Hz ((the sampling rate). Under 
these conditions the authors are correct in asserting that the 
aliasing test gives a positive result, which they believe to 
be incorrect. However, what the aliasing test actually in- 
dicates is that this signal is non-stationary. The particular 
interaction which the authors have constructed is not one 
for which the continuous-time selection criteria is met, i.e., 
the frequencies involved do not sum to zero. Samples of this 
signal do meet the discrete-time stationarity condition and 
so a non-zero bispectrum is possible in the outer triangle. 

One can look at these results in various ways. Our po- 
sition is that neither example constitutes a counter-example 
to the validity of the aliasing test in theory, though they both 
show that the test is limited in practice. The first example 
shows that there are aliased signals which the test does not 
see. This is obvious anyway since there are signals with 
zero bispectrum whose samples can be aliased. The sec- 
ond example shows that the term "aliasing test" must be 
restricted to stationary signals. As stated earlier, this re- 
striction is inherent in the definition of the bispectrum. 

6. CONCLUSIONS 

So, is this something for nothing? How can one get informa- 
tion about higher frequency amplitudes from what is usually 
thought of as Nyquist-limited data? The answer is of course 
that the assumption of stationarity is far from nothing. But, 



to exploit stationarity one must be able to perform the en- 
semble averaging indicated in the definition of the bispec- 
trum. This implies that one must either have an ergodic 
process or have access to sufficiently many sample paths. 

It is certainly possible that, in practice, the bispectrum 
can be usefully applied to signals for which there is no the- 
oretical justification. For such uses the aliasing test is silent. 
However, it is essential that a clear understanding of the 
fundamental properties of higher-order spectra be available. 
And the present authors believe that correct understanding 
of the outer triangle leads to deeper insight of the meaning 
of the bispectrum in general. 
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